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^ : ABSTRACT 

u 

-<^ ■ We present neutral hydrogen, ultraviolet, optical and near- infrared imaging, 

t^^ ■ and optical spectroscopy, of Minkowski's Object (MO), a star forming peculiar 

galaxy near NGC541. The observations strengthen evidence that star formation 
in MO was triggered by the radio jet from NGC541. Key new results are the 
discovery of a 4.9 x 10^ Mq double H I cloud straddling the radio jet downstream 
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•^ ■ from MO, where the jet changes direction and decoUimates; strong detections of 

■>!::j- ■ MO, also showing double structure, in UV and Ha; and numerous H II regions 
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and associated clusters in MO. In UV, MO resembles the radio-aligned, rest-frame 
O ■ UV morphologies in many high redshift radio galaxies (HzRGs), also thought to 

^ ■ be caused by jet-induced star formation. MO's stellar population is dominated 

by a 7.5Myr-old, 1.9 x 10^ M© instantaneous burst, with current star formation 
rate 0.52 MQyr^^ (concentrated upstream from where the H I column density is 
^ ■ high). This is unlike the jet-induced star formation in Centaurus A, where the 
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jet interacts with pre-existing cold gas; in MO the H I may have cooled out of 
a warmer, clumpy intergalactic or interstellar medium as a result of jet interac- 
tion, followed by collapse of the cooling clouds and subsequent star formation 
(consistent with numerical simulations). Since the radio source that triggered 
star formation in MO is much less luminous, and therefore more common, than 
powerful HzRGs, and because the environment around MO is not particularly 
special in terms of abundant dense, cold gas, jet-induced star formation in the 
early universe might be even more prevalent than previously thought. 

Subject headings: galaxies: jets — galaxies: starburst — stars: formation 



1. Introduction 

In a previous paper, twenty years ago, van Breugel et al. (1985) hypothesised that the 
starburst seen in the peculiar "Minkowski's Object" (MO; Minkowski 1958) in the cluster 
Abell 194 was triggered by the radio jet emerging from the nucleus of the nearby active 
galaxy NGC541. Only one other example of jet-induced star formation was known at the 
time, in Centaurus A (Graham & Price 1981; Brodie et al. 1983), and the phenomenon 
of jet-induced star formation was considered more of a curiosity than something worthy of 
general interest. In both cases the radio galaxies are at low redshift and of relatively low 
radio luminosity. Subsequent observations, with more sensitive detectors, larger telescopes, 
and better spatial resolution have shown that jet-induced star formation is more common 
than initially thought and examples have now been found in sources at low as well as high 
redshift, and with low (FR-I type) or high (FR-II type) radio luminosity (Fanaroff & Riley 
1974). 

For example, we now know that star formation in low redshift, low radio luminos- 
ity radio galaxies can occur in a variety of environments such as the central intergalactic 
medium (IGM) in X-ray clusters like Abell 1795 (McNamara 2002; O'Dea et al. 2004) and 
the interstellar medium (ISM) of their parent galaxies (e.g.. Wills et al. 2004). Some of 
the most dramatic jet-induced star formation, with star formation rates (SFR) as high as 
lOOOM0yr~^, is associated with very luminous radio galaxies at redshifts up to z ~ 4 (Dey 
et al. 1997; Bicknell et al. 2000). In all these cases the radio sources are embedded in rel- 
atively dense, and presumably clumpy, media as evidenced by X-ray halos, cold gas (H I), 
dust, bright line emission or a combination of these. Theoretical considerations suggest that 
the star formation may be induced by moderately slow (< 1000 kms^^) radiative shocks as 
they encounter IGM or ISM clouds (De Young 1989; Rees 1989). These qualitative models 
have been developed primarily for the powerful high redshift radio galaxies, where the star 
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formation is triggered by sideways expanding, radiative shocks in the relatively dense IGM 
/ ISM of the gaseous halos in proto-clusters or young parent galaxies (Bicknell et al. 2000; 
Klamer et al. 2004). 

However, little is known with certainty of the astrophysical conditions within the radio 
jets, including the composition of the jet fluid and its speed, as well as their environment, 
especially at high redshift. Better knowledge about this would be of great general interest, 
as it would help with modeling outflows from black hole accretion disks and the conditions of 
the ISM / IGM in forming galaxies and their proto-clusters. Numerical simulations are being 
used to explore the physical parameter range that is thought to be relevant to gain further 
insight. To help constrain such simulations MO is of special interest since in this case there 
is no obvious indication that the jet is propagating through some especially dense IGM or 
ISM environment. The Abell 194 cluster does not have a luminous "cooling" X-ray halo like 
Abell 1795, and there is no evidence for much cold gas outside the jet like in Centaurus A. 
Thus a better understanding of how the star formation in MO was triggered by the jet would 
also help in determining how much of a general phenomenon this might be. 

So why did the NGC 541 radio jet trigger star formation? Was it just a chance encounter 
with a gas rich galaxy that wandered into path of the radio jet, which was then "rejuve- 
nated"? This would have been a lucky shot, given that MO is the only peculiar system 
within the central ~ l,Mpc of Abell 194 and that its total (projected) diameter is the same 
as that of the radio jet. However, if there were an old stellar population underlying the 
starburst in MO, then this could not be ruled out. On the other hand, as the observations of 
Abell 1795 have shown, the interaction of jets with warm, clumpy gas in cooling flow clusters 
may well trigger the collapse of IGM clouds and subsequent star formation. 

The development of the new massively-parallel COSMOS numerical simulations package 
(Anninos, Fragile, & Murray 2003) at Livermore allowed Fragile et al. (2004) to simulate the 
collision of an FR-I type jet with a gas cloud in an effort to see whether star formation could 
be induced by the jet in the case of MO. They assumed that NGC 541 is surrounded by a 
multi-phase medium resembling galaxy cluster atmospheres (e.g. , Ferland et al. 2002); the 
simulations began with a hot "mother" cloud with a semimajor axis of lOkpc, a semiminor 
axis of 5kpc, density 0.1 cm~^, and temperature 10^ K. Within the mother cloud, denser, 
warm clouds were assumed to be embedded with typical sizes of 100 pc, densities of 10 cm~^, 
and temperatures of 10^ K. Studies of the proto-typical FR-I radio galaxy 3C31 (Laing & 
Bridle 2002) were used to estimate plausible parameters for the jet at the location of MO, 
15 kpc from NGC 541: a velocity of 9 x 10^ km s~^, a density of 10~^ cm~^, and a temperature 
of 10^ K. Fragile et al. (2004) concluded that it was indeed plausible that such star formation 
could be triggered by the driving of a radiative shock into such a cloud, causing the cloudlets 
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to compress, radiatively cool, and break up into numerous dense, cold fragments. These 
fragments survive for many dynamical timescales and are presumably precursors to star 
formation. 

Motivated by these results from numerical simulations, and by the availability of new 
observing facilities since the previous study, we find this an opportune time to revisit MO in 
an attempt to better understand the nature of jet / cloud interactions. We will first discuss 
the various new observations obtained, supplemented by archival observations. Then we 
proceed with an analysis of the spectroscopic and spectral energy distribution information, 
a discussion of the kinematics, morphology and time scales, and finally our conclusions and 
ideas for future work. We assume an fi^ = 0.27, Ha = 0.73, Hq = 71 km s~^ Mpc"^ 
cosmology (Spergel et al. 2003). At z = 0.019, MO and the Abell 194 cluster are both at a 
comoving radial distance of 78.7Mpc, and have a linear scale size of 0.374kpc arcsec"^. 



2. Observations 

2.1. VLA radio continuum data 

van Breugel et al. (1985) published 1.4 GHz Very Large Array (VLA) maps (taken in 
A-, B-, and C-array) of MO. These original radio data were re-reduced for the present 
investigation using standard routines in the AIPS software package. This was done in order 
to benefit from the cumulative improvements to the software over the years since 1985. This 
new radio continuum image is included in Figs. 1 and 2, and has a resolution of 3". We refer 
to the original paper for more details on the radio data. 



2.2. HST archival V and /-band WFPC2 images 

Although not the main target, MO was within the field of two short (460 s each) Hub- 
ble Space Telescope (HST) WFPC2/F555W and WFPC2/F814W observations of NGC541 
(Verdoes Kleijn et al. 1999), which are available in the archives. These are the highest spatial 
resolution images of MO, with a measured FWHM of 0.17" at F555W and 0.20" at F814W. 



-5- 



■40 - 



50 - 






.2 
C 

DJ 

Q 



1.0 - 



20 



30 - 



-0!*22W 




Or25'"49 



ie= 4?* 4e* 

Rlghl Ascension {J2QOO] 



45^ 



Fig. 1. — Coadded HST F555W / F814W image of MO, overlaid with radio continuum 
contours from the VLA. Contours are from to lOmJybeam^^ at intervals of 2mJybeam^^. 
The grey ellipse represents the 7'.'8 x 7'.'3 clean beam (PA = 36°) of the radio observations. 
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MO shows considerable structure with a total projected size which is approximately the 
same as the jet diameter (Fig. 1). The largest and brightest component (hereinafter referred 
to MO-North) is found along the northern edge of the jet. Its brightest portion, a 3kpc-long 
bar-like feature lies on the side closest to NGC 541 and extends at almost a right angle to 
the jet. A much fainter, linear feature (hereinafter referred to as MO-South) extends to the 
south across the jet. The central, high surface brightness portion of the radio jet bisects 
these two components just south of the position of the peak of the star formation as seen 
in the optical images (Fig. 1). Apart from the bar, MO-North contains several filamentary 
structures downstream in the jet, and is dotted with compact (mostly unresolved at the 
~ 75 pc resolution of HST) young star clusters and H II regions. 



2.3. Lick wide field i?-band imaging 

To investigate the local environment of MO, the central portion of the Abell 194 cluster 
was imaged by Michael Gregg on 2002 October 6, using the Prime Focus Camera (PFCam) 
on the Lick Observatory Shane 3-m telescope. Eight 600 s exposures were made using the 
Spinrad R filter, and five 600 s exposures were made in 5-band. The frames were reduced 
in the standard manner in IRAF, and coadded using offsets determined from individual 
astrometric solutions for each frame. The measured FWHM in the final image is 1.5" and 
the total field-of-view is 11' x 11' (Figs. 2 and 3). 
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Fig. 2. — Coadded PFCam {BR) image of the field of Abell 194 (color), overlaid with radio 
continuum (purple), H I data (dark blue), and Ha data (light blue). Minkowski's Object is 
clearly seen as the peculiar object at 1*^25™47'^, — 1°22'20", forming from the cloud of H I in 
the path of the jet from NGC541 (1^25™44^ -1°22'46"). The radio galaxy to the left of the 
image is 3C40 (1^26'^00^ -01°20'43"), associated with NGC547. The galaxy at 1^25™59^ 
-0r20'25", interacting with NGC547, is NGC545. 
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Fig. 3. — Zoom of figure 2 to show in more detail the region surrounding MO. 
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The image shows a large (~ 20 x 45 kpc) stellar bridge connecting the ellipticals NGC 541 
and the interacting galaxies NGC 545 / 547, which are associated with the large FR-I radio 
source 3C 40 (Simkin 1976). MO is located in this bridge, at least as seen in projection 
(Fig. 12), opening the possibility that the gas associated with MO has had its origin in some 
previous tidal interaction between these galaxies. The redshifts of these various systems 
are consistent: z = 0.0181 (NGC 541; Smith et al. 2000), z = 0.0189 (MO; this paper), 
z = 0.0178 (NGC 545; Smith et al. 2000), and z = 0.0182 (NGC 547; Trager et al. 2000). 



2.4. Narrow-band Ha imaging 

To determine the total extent of the star formation in MO and search for other possible 
emission-line features we obtained a narrow-band, dithered Ha image on 2004 December 
17 using "Imager" on the 2.3-m Advanced Technology Telescope (ATT) at Siding Spring 
Observatory. Three 200 s exposures of MO were made in the i?-band filter, and three 600 s 
exposures in a narrow-band filter centered at 668 nm (to match redshifted Ha in MO). Three 
30 s exposures in i?-band were made of the emission-line flux standard 205-26.7 (Dopita & 
Hua 1997), and three 100 s exposures in a narrow-band filter centered at 658 nm to provide 
the Ha flux calibration. The images in each filter were reduced in the standard manner in 
IRAF, using sky fiatfields appropriate to each filter. Once the individual images had been 
registered and co-added, an astrometric fit was performed for each coadd, and the broad- 
band images were registered and geometrically scaled so that the pixel positions of objects 
in the broad-band and narrow-band images matched to sub-arcsecond accuracy. 

Continuum subtraction was performed using the technique of Boker et al. (1999). Briefly, 
the reduced, registered broad-band and narrow-band images were rebinned to a size of 135 x 
135 pixels (each pixel representing 256 pixels of the original image). The intensity (in DN 
per second) of each pixel in the rebinned i?-band image was then plotted against that of 
the same pixel in the rebinned narrow-band image, excluding the brightest 500 pixels in the 
narrow-band image. A linear least-squares flt was performed, the slope of which represents 
the scaling factor between the two fllters. The unbinned i?-band image was then convolved 
with a Gaussian filter to match the PSF of the narrow-band image. It was then scaled by 
the slope of the fit obtained above, and subtracted from the unbinned narrow-band image. 

The scaling obtained for the standard star observations was 0.0438, while that for the 
MO observations was 0.0554. This is presumably due to differences in the narrow-band filter 
transmission curves (since the same i?-band filter was used for both fields), and the fact that 
the 668 filter samples a different part of the SEDs of the objects in the field to the 658 filter. 
The latter should be relatively unimportant due to the relatively small shift in wavelength 
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between the two filters, except in the case of strong emission hnes which may fall in one filter 
but not the other (e. g. , Ha in particular). Although this is the case for the main targets in 
each field, the good fit obtained when all binned pixels are included suggests that this effect 
is negligible, and that the main difference is due to the relative transmission of the filters. 
Therefore, we scale the Ha fluxes measured from MO by the ratio of the slopes (i. e. , the 
MO fluxes are multiplied by 0.791) to account for this. The resulting continuum-free Ha 
image is shown in Fig. 4 and has a resolution of 1.3" FWHM. 
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Fig. 4. — Left: Continuum-subtracted ATT Ha emission image of MO. Center: HST F555W 
(stellar continuum) image. Right: Keck NIRC K' image. All three images are aligned and 
show the same 1.0 arcmin^ region of sky. 
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The total flux integrated within an equivalent 16" diameter aperture (after correction for 
Galactic reddening - see section 3.1) is determined to be 8.6 x 10"^'* ergs s~^ cm~^, yielding 
an Ha luminosity of 6.6 x 10^° ergs s~^. This is approximately twice the value obtained 
by van Breugel et al. (1985) - some 14% of this increase is attributable to the different 
cosmology used in their paper, and the remainder is presumably due to improvements in 
both the sensitivity and calibration of the current Ha observations. 

The Ha image is similar in its morphological features to the optical continuum images, 
with a strong concentration of H II regions and their associated clusters in the linear bar-like 
feature facing the incoming jet, and more diffuse filaments of emission on the downstream 
side. These are more conspicuous in Ha than in the continuum, suggesting that the clusters 
in this region are somewhat younger than in the bar. MO-South is more clearly defined in 
Ha. The overall morphology is highly suggestive of a jet-cloud interaction, with subsequent 
star formation triggered by the jet. 



2.5. Keck ESI i?-band and NIRC J and i^'-band imaging 

To complement the HST V and /-band observations, and to provide spectral energy 
distribution data over the widest possible range, we obtained Keck i?, J and i^'-band images. 

In the i?-band, a 75 s image centered on MO was obtained with the Echellette Spectro- 
graph and Imager (ESI; Sheinis et al. 2002) on Keck II on 2002 September 7, and calibrated 
and reduced in IRAF in the standard manner. The FWHM of the seeing was 0.7" and the im- 
age is shown in Fig. 5. This image clearly shows the existence of a fainter structure extending 
about 30" downstream of the bright bar of MO-North, which was also faintly detected in the 
Ha image. 
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Fig. 5. — Keck ESI i?-band image of MO (top) with elliptical galaxy models subtracted 
(bottom) as described in the text. Both images are displayed using the same logarithmic 
scaling. The arrows in the bottom panel illustrate the filamentary feature downstream of 
the central bar, discussed in section 2.5. They are not intended to suggest that this feature 
necessarily formed due to flow of material from the central regions of MO. 
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600 s exposures in J and K' were obtained with the Near Infrared Camera (NIRC) 
on Keck I on 2004 December 17. A 5-position dither, with dither spacing 5" and 60s per 
exposure (3 x 20 s in J, 4 x 15 s in K'), was used. Between each on-source dither position, a 
40" nod to a nearby blank area of sky was made, and an additional 60 s exposure off-source 
was made (yielding the same dither pattern as for the on-source exposures). The on-source 
images were registered and combined. The off-source images were combined with object 
masking and without registration to form a sky image, which was scaled and subtracted 
from the on-source image. The final reduced images cover a field of 100" x 100" with a 
FWHM of 0.5". 

The i^'-band image is shown as the third panel in Fig. 4. The K'-band morphology 
very closely resembles the optical F555W image taken with the HST. This shows that there 
cannot be a large old stellar population present in MO, and that the i^'-band light is probably 
dominated by young supergiants. 



2.6. GALEX imaging 

The field containing the NGC541 / 545 / 547 galaxies in Abell 194 was observed as 
part of the GALEX (Martin et al. 2005) Nearby Galaxy Survey (Gil de Paz et al. 2006) 
in the near-UV (NUV, 1750 - 2800 A) and far-UV (FUV, 1350 - 1750 A) passbands on 
2003 October 23. Exposure times were 1700 s in each band. The data were reduced and 
photometrically calibrated using the GALEX pipeline (Morrissey et al. 2005). The 5(7 AB 
limiting magnitude in both bands is tti ~ 23 corresponding to a surface brightness sensitivity 
of ~26.5 magarcsec"^ averaged over the PSF. The FWHM in the NUV and FUV images is 
5.6" and 4.0" respectively. 

MO was prominently detected in both bands and shows similar overall structure to the 
optical. Ha and near-IR images, even given the relatively low resolution of GALEX. The 
GALEX field of view is large and covers the central 1.2° of the Abell 194 cluster, yet MO is 
among the brightest resolved UV structures in this entire field. 

The images were overlaid with the VLA radio continuum and H I observations (Sec- 
tion 2.7), as shown in Fig. 6. In both UV bands MO appears as a double source, with a clear 
separation of MO-North and MO-South straddling the central axis of the radio jet, with MO- 
North being the brightest component. Both components are clearly extended downstream, 
paralleling the radio jet, as in the optical and near-IR images. 
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Fig. 6. — Left: GALEX NUV image, overlaid with H I contours from the VLA. The NUV 
image is shown with a logarithmic scaling to emphasise faint structures; note that the blob 
at 0l'^25™50'^ is a foreground star. H I contours are from 0.5 - 6.5mJybeam~^ kms~^, at 
intervals of l.OmJybeam"^ kms^^. The grey ellipse represents the 17'.'9 x 15'.'4 clean beam 
(PA = 30°) of the H I observations. Right: GALEX FUV image (shown with a logarithmic 
scaling), overlaid with radio continuum contours from the VLA. Radio contours are from 2.0 
- 17.0mJybeam~^ at intervals of 3.0mJybeam~^. The grey ellipse represents the 7'.'8 x 7'.'3 
CLEAN beam (PA = 36°) of the radio observations. Note that the two panels are not shown 
to the same scale; the relationship between the radio jet and the H I is better seen in Fig. 8. 
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The total observed near-UV flux from MO was determined to be 1.78 x lO^^^ergs s~^ 
cm~^ A~^, and tlie far-UV flux was 3.63 x 10"^^ ergs s~^ cm~^ A~^, yielding values corrected 
for Galactic extinction of 2.57 x lO^-'^^ergs s~'^ cm~^ A^^ and 5.11 x 10^"*^^ ergs s~^ cm~^ 
A^^ respectively. This is in good agreement with a flux of 2.96 x 10^^^ ergs s^^ cm^^ A~^ 
at 1350 A from an archival International Ultraviolet Explorer spectrum (K. Chambers, priv. 
comm.). We also extracted photometry for three smaller, rectangular apertures, centered on 
MO-North, the NE fllament downstream from MO-North, and MO-South (Fig. 7). These 
GALEX photometry data allow us to put strict limits on the stellar age distribution of MO, 
as we shall discuss in Section 4.2. 
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Fig. 7. — Binned ESI-R image of MO, overlaid with boxes showing the regions used for the 
three-aperture photometry (see text). 



-18- 



2.7. The Radio H I VLA observations 

To search for cold gas associated with MO and its surroundings, we performed H I obser- 
vations with the VLA in its C- array configuration on 2002 April 2. The total integration time 
on source was 4.3 h. A 3 MHz bandwidth was used, using 128 spectral channels, covering the 
heliocentric velocity range 5305 kms^^ to 5815 kms~^(the full observing band is 600 kms~^, 
but on either side of this band are a few channels that are not usable). The data were cal- 
ibrated following standard procedures using the MIRIAD software package. A spectral-line 
data cube was made using robust weighting (robustness = 0.5) giving a spatial resolution 
of 18" X 15" (PA = 29°). Banning smoothing was applied to the datacube, resulting in a 
velocity resolution of 10.3 kms~^. The noise in the final datacube is 0.54mJybeam~^. 

Within the entire primary beam (0.5°) the only H I that we detect is from MO (Fig. 8), 
but this may be due to our relatively narrow band covering only ~ 500 kms~^ (the velocity 
range of the cluster is surely larger than that). Indeed, a short test observation with the 
WSRT in August 2004 of MO with a wider band of 10 MHz covering 2000 km s"^ (using 1024 
channels) showed at least 3 other H I objects, apart from MO, at velocities just outside the 
band used with the VLA for the MO observations. A more detailed observation of Abell 194 
with the WSRT would likely reveal more objects. 
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Fig. 8. — Map of the H I cloud adjoining MO, with radio continuum contours from the VLA 
overlaid. The H I cloud is shown with a linear scaling from 0.3 - TmJybeam^^ kms^^, corre- 
sponding to the greyscale wedge at right. Continuum contours are from 1 - 17mJybeam~^ 
at 2mJybeam~^ intervals. The H I beam is shown as a grey ellipse (lower left), and the 
continuum beam is shown as a white ellipse (upper right). 
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The H I emission is dominated by two "blobs" sitting on either side of the jet and is 
offset from MO, with the brightest H I ~ 10" (3.6 kpc) downstream from the peak of the Ha 
emission (the main star formation region). The position of these maxima is coincident with 
the downstream extensions of MO-North and MO-South noted in the GALEX image. In 
addition, two faint H I filaments are found extended along the edges of the jet. One of these 
lies to the north of MO-North - upstream in the sense of the jet motion. The other forms 
a continuation of MO-South downstream. The jet itself changes direction and decollimates 
downstream from the H I clouds, which suggests that it is the jet-gas interaction which has 
disrupted the jet. 

The H I kinematics, along an east-west slice (see Fig. 9) show that most of the H I 
gas is relatively quiescent with a average velocity vhi ~ 5640 kms~^, or z = 0.0188. This 
is in excellent agreement with the redshift z = 0.0189 of the star-forming region obtained 
from the LRIS spectroscopy. To make this plot, the data were integrated in declination over 
about 35". This was done to show the faint low- and high- velocity wings. There is a smooth 
north-south velocity gradient and the North and South H I blobs have a velocity difference 
of about ~ 50 kms^^. 
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Fig. 9. — Velocity slice through the MO H I map, showing velocity shear near the outer 
boundaries of the H I cloud and jet. Contour levels are 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 mJybeani"^, 
and encode the same information as the grey scale. 
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We derive a total H I mass for MO of 4.9 x 10* Mq. This was done assuming the H I 
is optically thin. Then the mass follows from Mm = 2.35 x lO^FD"^, where F is the flux 
integral in Jy kms~^ (i.e. the integral J S dv where S is the detected flux in a channel and 
V is velocity) and D is the distance to MO in Mpc. 

The mass of H I is approximately evenly divided between the two main blobs straddling 
the jet (Fig. 8). The H I column density through the center of each of the blobs is hhi ~ 
4.3 X 10^*^ cm"^ and Uhi ~ 4.7 x 10^*^ cm"^ for the northern and southern blobs respectively. 
The total (projected) extent of the H I emission is approximately 275 kpc^. 



2.8. Keck LRIS spectroscopy 

Any potential spatial variations in the stellar populations or ages are best studied 
through spatially-resolved spectroscopy. To investigate any spatial variations in ionization 
and anomolous kinematics across MO we obtained a high signal-to-noise longslit spectrum 
using the Low Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) on Keck I on 2004 
January 20. The D560 dichroic was used with a slit width of 1.5". The total exposure time 
was 1800 s. On the blue side, a 300 line mm~^ grism, blazed at 5000 A was employed, giv- 
ing 1.43 Apixel"^ and spectral resolution 10. 7 A. On the red side, a 400 line mm~^ grating, 
blazed at 8500 A was used, giving 1.86 Apixel"^ and spectral resolution 6.8 A. The data were 
reduced in IRAF, with cosmic ray removal, sky line subtraction, wavelength calibration and 
flux calibration performed in the standard manner. The blue- and red-side spectra were 
combined using scombine in IRAF, and the resultant pixel scale was 1.41 Apixel"^, so the 
resulting spectra were smoothed with a 3 pixel boxcar fllter to improve signal-to-noise. Spec- 
tra were extracted using several apertures and were corrected for Galactic extinction using 
E{B -V) = 0.045 and Ry = 3.1 (Schlegel et al. 1998). 
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Fig. 10. — The LRIS slit position (left), overlaid on the F555W image of MO (the width of 
the slit is accurately shown). Numbered sub-apertures are the same as shown on the section 
of the 2-D spectrum (showing H/3 and the [O III] doublet) at right, and the fluxes in table 1. 
The "full slit" spectrum was extracted in an aperture covering subapertures 1-8. 
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A large, 21"-long aperture spectrum (just covering subapertures 1 - 8 as shown in 
figure 10) was extracted to obtain an integrated spectrum encompassing all of the emission 
from MO along the 1.5"-wide slit. The extracted spectrum is shown in figure 11. From this, 
the redshift of MO was determined to be z = 0.0189 (similar to the value found by van 
Breugel et al. 1985). The spectra were Doppler corrected to zero redshift. The measured 
emission line properties integrated along the slit, and for various sub- apertures described 
below are presented in Table 1. 

Fluxes were measured from Gaussian fits to the lines, with 1-sigma errors estimated from 
Monte Carlo noise simulations (noise estimates and gain were determined from measurements 
of the spectrum). 

The strengths of the Balmer absorption lines (Table 2) were measured by fitting a 
Gaussian to the wings of the absorption underlying the Balmer emission components. Quoted 
fluxes are the amount of absorbed flux based on a fit to the continuum after subtraction of a 
Gaussian model for the emission component. The emission line fluxes in Table 1 are corrected 
for underlying absorption in the line cores. 



Table 1. Equivalent widths and Galactic-reddening-corrected line fluxes from LRIS spectroscopy 

Line Xrest (A) EWo (A) Flux (H/3 = 100): 

Full slit Full slit Ap 1 Ap 2 Ap 3 Ap 4 Ap 5 Ap 6 Ap 7 Ap 8 



on 


3727 


63.6 ±0.4 


359.0 


221 


406 


449 


305 


454 


362 


289 


240 




[Ne TIT] 


3869 










40 


28 




31 








H8 


3889 


1.0 ±0.3 


7.9 




















m 


4102 


1.1 ±0.3 


9.7 




















H7 


4340 


4.7 ±0.3 


40.8 


51 


31 


33 


30 


31 


30 








H/3 


4861 


12.4 ±0.3 


100.0 


100 


100 


100 


100 


100 


100 


100 


100 




III 


4959 


6.7 ±0.3 


60.7 


115 


115 


59 


85 


49 


36 


67 


87 




[0 TIT] 


5007 


19.4 ±0.3 


172.4 


333 


333 


180 


252 


169 


114 


211 


200 




He I 


5876 


1.7 ±0.4 


9.4 




11 


17 


16 


11 


7 






to 


[01] 


6300 


2.5 ±0.3 


14.1 




22 


28 


20 


11 


16 




46 




01 


6363 


1.0 ±0.3 


5.2 




7 




7 












Nil 


6548 


2.9 ±0.3 


13.6 




15 


7 


23 


16 


10 








Ha 


6563 


58.6 ±0.3 


278.7 


226 


306 


335 


396 


318 


238 


350 


390 




Nil 


6583 


8.3 ±0.3 


39.5 


14 


44 


54 


59 


46 


36 








[SII 


6716 


17.1 ±0.3 


78.7 


31 


94 


105 


114 


81 


74 


100 


120 




SII 


6731 


12.2 ±0.3 


55.0 


18 


55 


75 


73 


57 


53 


50 






-^H/3 






789 ±2 


84 ±3 


124 ±3 


57 ±3 


116 ±3 


140 ±3 


134 ±3 


18 ±3 


15 ±4 





Note. — H/9 line fluxes are given in 10 ergs s cm 



to 
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Table 2. Absorption lines from LRIS spectroscopy, corrected for Galactic reddening 
Line Xrest (A) Flux^ EWq (A) 



HIO 


3798 


218 ±2 


3±1 


H9 


3835 


262 ±2 


3±1 


H8 


3889 


223 ±2 


3±1 


m 


4102 


260 ±2 


3±1 


H7 


4340 


363 ±2 


4±1 


H/3 


4861 


247 ±1 


3±1 



Note. — ^Line fluxes are given in 10 ^^ ergs s ^ cm 
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Additionally, spectra were extracted in eight subapertures, chosen to coincide with re- 
gions of the 2-D spectrum and images which appeared to have relatively homogeneous prop- 
erties across them (see Fig. 10). The properties of the measured emission lines are also shown 
in Table 1. There was insufficient signal-to- noise in the absorption lines in the subapertures 
to accurately determine absorption line strengths at these positions. 
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Fig. 11. — The combined blue and red side LRIS spectra (1800 s), smoothed with a 3-pixel 
boxcar filter. 
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Line ratios which depend upon the measurement of hues on both the red and blue sides 
are less reliable for the subapertures than for the integrated spectrophotometry, due to the 
difficulty of measuring exactly the same region of the 2-D spectrum in apertures of only a 
few pixel length. This is presumably part of the source of the variation in the estimated 
E{B — V) seen in Table 3, below. Line ratios determined using lines from only the red or 
blue sides are unaffected by this additional source of uncertainty. 

Since the slit is oriented north - south, the observations will to some extent be affected 
by atmospheric differential refraction (Filippenko 1982). This was minimised by observing at 
relatively low airmass {Z < 1.20), but dispersion still has the effect of offsetting the physical 
regions of the source where lines at the red and blue ends of the spectrum are measured 
(as badly as 0.9" in the case of [O II] and [S II]). In the case of ratios such as Ha /H/5, the 
differential refraction is at worst around 0.36", i.e. , around 25% of the width of the slit. 
Inspection of figure 15 suggests that the extent of the extended emission compared to the 
slit width is such that a shift of this magnitude ought not to produce wild fiuctuations in 
line ratios, especially since the seeing at the time of the LRIS observations (0.84", some 5 
times worse than that of the F555W image in figure 15) will tend to wash out features at 
smaller scales. 

In an extreme case where the fiux of Ha or H/5 was different by a factor of two in the 
25 % of the slit discussed above, this could lead to errors in the computed extinction values 
of E{B — V)s ~ 0.08 (see section 3.1). The fact that the calculated value of E{B — V)s for 
the whole slit is so close to zero (formally, E{B — V)s = —0.01) rather than, for example, 
strongly negative, suggests that these assertions are reasonable. This effect is presumably 
somewhat more important on smaller scales (e. g. , the individual slitlets) and so may be an 
additional source of uncertainty (along with the difficulty of measuring the same region on 
the red and blue sides of the spectrum, discussed above) in the E{B — V) values in table 3. 
Since in practice, for the whole slit, we do not expect variations in ffux of a factor two on 
scales smaller than the atmospheric seeing, we assume that atmospheric dispersion is not a 
major source of error in our analysis. 

The astute reader will note that the total H/5 ffux quoted for the "full slit" in table 1 
is around 15% larger than the sum of the ffuxes from the individual subapertures. This is 
again due to the difficulty of placing individual subapertures to sub-pixel precision, but as 
noted above does not affect ratios as measured in individual apertures. 

The line ratios observed in these spectra confirm that the line emission is due to H II 
regions throughout MO (Veilleux & Osterbrock 1987; Dopita et al. 2000). In particular, 
there is no evidence of shock-excited gas, which would reveal itself by very large [N II] / H/5 
or [O II] / B.(3 ratios (Dopita et al. 1984). Nor does the [O I] /Ha ratio become sufficiently 
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large to put it in the shock-excited classification (Dopita & Sutherland 1995). Furthermore, 
the kinematics show that the ionized gas is quiescent. There is no evidence for large velocity 
gradients, at least across the jet in the north-south slit direction. However, since shock- 
excited gas is more likely to reveal itself at low fiux levels, it might be of interest to obtain 
high signal-to-noise spectroscopy along the jet direction. 



3. Derived Quantities 

3.1. Reddening 

The measured Balmer line ratio Ha / B.f3 from the full LRIS spectrum (after correction 
for Galactic reddening), was found to be 2.790 ± 0.008. Following Calzetti et al. (1994) and 
Calzetti (2001), this is consistent with an intrinsic (continuum) reddening oi E{B — V)s ~ 0. 

We also calculated values oi E{B — V)s for the subapertures (see table 3). These show 
some variation, including values as negative as —0.10. We believe these determinations 
to be less reliable than the full spectrum value for reasons discussed in section 2.8; errors 
introduced by atmospheric dispersion and by not measuring the same part of the spectrum 
in the red and blue sides can probably account for any deviations from E[B — V)s = in 
the subapertures. 

We likewise used the (Galactic reddening corrected) broad-band UV fiuxes from GALEX 
to calculate f3, the UV continuum slope (Calzetti et al. 1994), and obtained a value of 
f3 = —1.734, again yielding E{B — V)s ~ 0, in excellent agreement with the line ratio 
results. 

These more accurate measurements imply that there is essentially no dust obscuration 
along the line of sight to the emission line regions of MO and supersede those of van Breugel 
et al. (1985) who obtained E{B — V) = 0.27. It is possible, however, that some of the 
variation in calculated values oi E{B — V) in the subapertures is in fact due to the presence 
of localised dust. If the dust is mixed in with the emission line gas and stars to some extent 
(rather than a foreground screen, which is the assumption in the extinction calculations) 
then dust mass can be somewhat higher for a given E{B — V) than otherwise. 



3.2. Metallicity 

As remarked above, the emission-line ratios are typical for those of H II regions, show- 
ing that the emission is dominated by star forming regions and so confirming the earlier 
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conclusions by van Breugel et al. (1985). The high quahty of our data allows us to obtain 
additional, more accurate diagnostics. 

We determined the metallicity of MO using the method outlined in Kewley & Dopita 
(2002). This yields a value of log{0/H) + 12 = 8.6 ± 0.1, equivalent to a metallicity Z ~ 
0.5 Zq. The mean ionisation parameter was determined to be g ~ 3 x 10^ cms^^ {U = 10^^). 
This is within the normal values expected for extragalactic H II regions (Dopita et al. 2000). 



3.3. Star Formation Rates 

Applying the relation between Ha luminosity and star formation rate (SFR) (Dopita & 
Ryder 1994; Kennicutt 1998; Panuzzo et al. 2003), 



(7.0-7.9) X 10" 



Mq jt~^ J yergs"^ 

we used the measured Ha flux from the ATT images to obtain a total SFR within MO 
of ~ 0.5 Mq yr~^. This estimate is sensitive to the form of the initial mass function assumed 
for high mass stars, and therefore has a uncertainty of order 30%. 

The local star formation rate increases, as expected, in the subapertures near the central 
blobs, where it is up to an order of magnitude higher than in the subapertures in the outer 
regions of MO . The measured star formation rates are given in Table 3, assuming the globally 
determined result of no intrinsic dust. 
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Table 3. Derived quantities from the observed emission lines from the LRIS subapertures 



Quantity 




Ap 1 


Ap2 


Ap3 


Ap4 


Ap5 


Ap6 


Ap 7 


Ap 8 


E{B-V) 




-0.10 


0.03 


0.06 


0.13 


0.04 


-0.08 


0.08 


0.13 


Subaperture width (") 




1.89 


5.13 


1.35 


1.22 


1.35 


2.97 


1.78 


2.43 


Subaperture size (kpc^) 




0.397 


1.076 


0.283 


0.256 


0.283 


0.623 


0.373 


0.510 


SFR(Ha) 10-3 Moyr-i 


kpC-2 


2.90 


2.14 


4.09 


10.90 


9.56 


3.12 


1.02 


0.70 


;S II] 6716 / 6731 




1.73 


1.71 


1.40 


1.56 


1.43 


1.39 


2.00 




logW^ 




-2.6 


-3.0 


-3.2 


-3.0 


-3.0 


-3.0 







^Determined from the [N II] / Ha diagnostic (Kewley & Dopita 2002), assuming constant 
metalhcity oi Z = 0.5Zq 
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3.4. Electron Densities 

The density-sensitive [S II] A6716/A6731 doublet ratio was measured to be 1.431 ± 0.01 
in the full spectrum. No other electron density diagnostics were able to be determined from 
our LRIS spectra, due to insufficient spectroscopic resolution and / or sensitivity. Using the 
five-level atom model of de Robertis et al. (1987), and assuming an electron temperature, 
Te = 10^ K, this yields an electron density of order 1 - 10cm~^, averaged along the slit. 
However, for electron densities lower than approximately 100 cm~^, corresponding to values 
of this ratio > 1.32, this line ratio becomes rather insensitive to electron density, and our 
measured electron density is therefore highly uncertain. Measurements of the ratio in the 
individual subapertures yield values ranging from 1.39 - 2.00 (Table 3). Given that the 
errors on the ratios are rather high, due to low signal-to-noise in some subapertures, and that 
the S II lines are almost coincident with strong telluric lines whose subtraction introduces 
additional uncertainties, all we can conclude with certainty is that the electron densities are 
less than 100 cm~^. 



3.5. Ages of the Clusters 

As Gonzalez Delgado et al. (1999) have shown, the strengths of the He I lines relative 
to H/3 are sensitive to the age of an H II region. We created our own photoionization models 
(appropriate to the metallicity of MO determined above) using mappings lllr (Sutherland 
& Dopita 1993; Dopita et al. 2002) in combination with Starhurst99 (Leitherer et al. 1999) 
stellar synthesis instantaneous burst models, with a Salpeter IMF with lower mass limit 
of IMq and upper mass limit of 100 M©. The use of a Salpeter IMF is appropriate for 
the modelling of both the ionizing flux and stellar continuum for a young (< 10^ yr) recent 
starburst. The neglect of the < 1 Mq stars makes no difference to the results, as they 
contribute negligible luminosity at these early times. We flnd that the He I A5876 / H/9 
line ratio remains near 0.1-0.12 until 4Myr, after which it declines steeply, becoming < 0.05 
after 5 Myr. The measured line ratios in table 1 are therefore consistent with H II region 
ages of < 5 Myr, for all sub-apertures. 

The [O III]A5007 / H/? ratio is even more sensitive to age. In our models, it starts at 
around 4.5, declines to ~ 1.2 by 2.5 Myr, to 0.5 by 3.0 Myr, and is negligible thereafter. We 
can therefore conclude that all the observed H II regions are younger than ~ 3 Myr. 

Comparison of the equivalent widths (EWs) of the Balmer nebular emission lines (ta- 
ble 1) with flgure 16 of Gonzalez Delgado et al. (1999) rules out continuous starburst models. 
The measured H/3 and H8 EWs both yield an instantaneous burst age of 6 Myr, and H(5 yields 
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an instantaneous burst age of 8 Myr. Comparison of the measured EWs with figures 3 and 4 
of Gonzalez Delgado et ah (1999) shows that ages greater than around 10 Myr are excluded. 

The discrepancies between the ages obtained from [O III]A5007 / H/3 and those from 
the EWs (discussed above) and SED fitting (Section 3.8) may be due to a combination of 
some ongoing star formation, or possibly some contribution to the [O III] luminosity from 
shocked gas. Spectra with higher spatial resolution and signal-to-noise might help resolve 
this. 



3.6. Number of Ionizing Photons 

Taking the Galactic reddening corrected ratio Ha / H/? = 2.79 determined from the 
full-aperture spectroscopy, and the Ha fiux for the whole object from the imaging, we derive 
an H/9 luminosity of 2.37 x 10*^*^ ergs s^^. Comparing to Gonzalez Delgado et al. (1999) table 
7, and assuming an age of 6 Myr (from the emission line EWs), we obtain a stellar mass of 
1.9 X 10^ Mq, and log Q* (the number of Lyman continuum photons) of 52.7. If we were 
instead to use the age of 3 Myr obtained from the H II region spectroscopy, this would yield 
a stellar mass of 1.5 x 10^ Mq, so this determination is quite sensitive to age. The value of 
log Q* depends only on the H/5 luminosity and so is independent of age. 

From Stasinska & Leitherer (1996), 

W = 2.8 X 10-2°(10VTe)2/3(g,ne2)V3 

and with the measured values U = 0.001, n = 10 cm"'^, log Q = 52.7, and assuming 
Te = 10'^ K, we obtain a nebular gas filling factor e ~ 0.01. 

3.7. Photometric Properties 

The NUV, FUV, ESI-R, F814W, J and K' images were registered to the F555W im- 
age using the individual astrometric solutions. The F555W, F814W, J and K' images were 
convolved to match the 0.72" seeing of the ESI-R image. We refrained from convolving all 
images to match the 5.6" FWHM of the NUV image since this would result in such a sig- 
nificant loss of spatial resolution for the majority of the images. Photometry was performed 
in an 8" radius circular aperture to enclose all the emission from MO. Since MO is resolved 
in the GALEX images (Section 2.6) we also obtained photometry for three sub-apertures 
(Fig. 7) centered on MO-North, the NE filament, and MO-South (Section 4.2). 
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The elliptical galaxies in the ESI image (Fig. 5) were modelled using the tasks ellipse 
and BMODEL in IRAF, and the best fitting models were subtracted from the original image. 
NGC 541 is difficult to fit in this image because it extends beyond the boundary of the frame 
to the right, but we were able to do a reasonable job of subtraction by carefully choosing 
and refining the fit parameters (see Fig. 5). 

After subtraction of the elliptical galaxies, the bridge of stars joining NGC 541 and 
NGC 545 / 547 (see also Fig. 12) can be more clearly seen and more accurately measured. 
The bridge is around 24.5 mag arcsec"^, and extends over around 115"x 45" (~ 700 kpc^); 
its total R-band luminosity is around 3 x 10^ Lq. From Nikogossyan et al. (1999), the gas 
mass within a 0.8 Mpc radius of the center of the cluster is 1 x 10^^ Mq, so we expect (by 
simple scaling, in the lack of more detailed information about the gas distribution) the gas 
mass in the bridge to be approximately 3 x 10^ Mq. 
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Fig. 12. — PFCam i?-band coadd, cropped to highlight the faint bridge of stars which 
connects NGC541 with the interacting pair NGC545 / 547. MO hes in this bridge, close to 
NGC541. 
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The bridge is approximately 18.7 mag in the aperture used for the broadband photom- 
etry (i. e. , approximately 15% as bright as MO) and the portion of the NGC 541 model 
enclosed by this aperture is 55I20 % ^^ bright as MO. Hence about 40% of the light in 
the aperture comes from NGC 541 and the bridge. Since a circular aperture was used, this 
includes parts of the aperture where MO is relatively faint, and parts where it is relatively 
bright, and so cannot be construed to be typical for the spectroscopic apertures or the in- 
dividual pixels of the image. However, we consider the contribution from the bridge and 
elliptical when performing spectral synthesis modeling. 



3.8. Spectral Energy Distributions 

We computed SEDs using the seven bands of integrated (8" radius aperture) Galactic 
reddening corrected photometry for MO, converting magnitudes to fluxes in ergs s~^ cm~^ 
A^^ using published zeropoints, and normalising the fluxes to the flux measured in the 
F555W fllter (approximately V"-band) to facilitate model comparisons. A set of model spectra 
were created using Starburst99, with parameters as described in section 3.5 . These models 
were redshifted to match the redshift of MO, and convolved with the fllter response functions, 
resulting in a set of fluxes for each model. These were normalised to the F555W flux, and the 
best match between these fluxes and those from the observed photometry were determined 
by a chi-squared minimisation procedure. 

In flgure 13, we plot xl = X^/'^j where z/ = 6 is the number of degrees of freedom (since 
seven fllters are used, and the normalisation takes away one degree of freedom). The best 
fltting (single-age) model was a 7.5 Myr old SED with xl = 0.22. Models older than 14Myr 
flt the photometry at V through K' quite well, resulting in a broad minimum in the x^ fit 
(2 ^ X^ ^ 4). However, these older models do a poor job of reproducing the UV flux, and 
hence have higher x^ than the 7.5 Myr model, which results in a sharp minimum in x^, and 
reproduces the UV fluxes much better. A Ax^ analysis, following Avni (1976) (with the 
number of "interesting parameters" equal to one, i.e. , the template age), yields formal 68% 
confidence limits of 7.47 ± 0.12 Myr (interpolating xl values between the 0.25 Myr template 
age spacing). The secondary minimum at 14.02 Myr {xl = 0.80) is statistically a worse fit 
at > 90% confidence. The associated x^ probability for the primary minimum is 0.97. 

The two disjoint regions of low chi-squared are due to the similarity of the SEDs at 
7.5 and 14.0 Myr. Between these ages, the UV and infrared fiuxes increase above the model 
values due to the dominance of red supergiants in the models. As mentioned above, older 
models (> 25 Myr) fit the optical and infrared points reasonably well, but do not do a good 
job of fitting the UV points. 
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In an attempt to account for an underlying (presumably old) population from the bridge 
and NGC541 (see section 3.7) we ran the x^ minimisation procedure again, this time adding 
in a 1 Gyr population at the 40 % level (the number of degrees of freedom is unchanged 
since the normalisation of this component is fixed by the photometry from section 3.7). The 
best fit for the remaining component was still the 7.5 Myr model, with x^ = 1.02 (slightly 
higher than for the single population model, but with a x^ probability 0.41). This yields 
68 % confidence limits of 7.501q'}8 Myr, with a secondary minimum at 13.94 Myr (x^ = 1.23) 
statistically distinct from the primary solution only at > 68 % confidence. 

The old population fraction was increased in steps of 1 %, and the x^ procedure run 
over again. The x^ minimum continued to increase (as the UV data points were less well fit), 
but the 7.5 Myr model was still selected. At an old population fraction above 50 %, the best- 
fitting young population begins to fluctuate between models and the fit becomes increasingly 
poor. It seems reasonable, therefore, to conclude that less than 20 % (and possibly as low as 
%) of the light (although possibly a larger percentage of the stellar mass) in the aperture 
is due to an old population in MO itself, as suggested by the photometry in section 3.7. 

We used a technique similar to that employed by Johnston et al. (2005) to obtain an age 
map for the galaxy. The registered, seeing-convolved images were binned 4 by 4 pixels. The 
model fitting routine, as described above, was run on the SEDs extracted from individual 
corresponding binned pixels in each band, producing a best-fitting age for each pixel (again 
assuming a single- age instantaneous starburst). The resulting pixel age map is shown in 
Fig. 14. 
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Fig. 13. — Left: A selection of synthetic spectrum models (solid lines) and the photometry 
computed by convolving those models with the filter transmission curves (colored open cir- 
cles), is compared to the observed (Galactic reddening corrected) photometry (black points). 
The vertical errorbars on the black points represent measured photometric errors, and the 
horizontal bars represent the filter effective bandwidths. The best fitting (7.5 My) model is 
shown in red, along with a badly fitting model (9.5 My in green) and another model which 
provides quite a good fit, i. e. low chi-squared (14.0 My in blue). Right: Reduced chi-squared 
for the model fits, as a function of time since the instantaneous burst in the spectral tem- 
plate (solid line). There is a clear minimum around 7.5 My, and another broad, but less 
deep, minimum for ages older than 14 My. The old models fit the photometry at V through 
K' quite well, but do a poor job of reproducing the UV fiux. Adding a 1 Gyr-old population 
at the 40 % level produces x^ values shown as the dotted curve. 
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Fig. 14. — Left: Pixel age map for MO. A range in ages is seen in MO, color coded according 
to the wedge at right. Dark blue pixels were rejected due to their large x^ (meaning they 
were poorly fit). Red pixels were best fit by models of > 9 Myr (this includes ages up to and 
including 1 Gyr). White contours are from the binned HST F555W image, and are overlaid 
to show the relationship between the morphology of MO and the clustering of pixels of the 
same age. Low signal to noise results in some spurious pixels, but broadly it appears that 
pixels of the same age cluster together, with the oldest populations in the central bar of MO, 
and younger populations in the outer regions. This picture is supported by binning pixels 
in these regions to improve signal-to- noise (see text). Right: Histogram of the best fitting 
pixel ages. Most of the older pixels are in the nearby old elliptical NGC541. The peak at 
40 My includes the 1 Gyr old pixels. 
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Reassuringly, pixels of the same age cluster together. There appears to be a gradient in 
MO, from pixels with ages around 8.5 Myr at the center, to pixels with ages around 7.0 Myr 
at the edges. The edge of NGC 541 at the lower right of the frame, and also the nearby small 
background galaxy, are fit with a 1 Gyr old population (the oldest population considered in 
our fits). We are confident, therefore, that the pixel ages are representative of the dominant 
population. Since young stellar populations tend to dominate the integrated light of a galaxy, 
we cannot rule out an underlying old population, and indeed a small number of pixels in MO 
itself are best fit with older populations. However, this may be simply due to low signal-to- 
noise in the UV data in these pixels, which provides less discrimination between the primary 
and secondary minima in the x^-age plot (the primary minimum tends to remain relatively 
sharply defined, whereas the secondary minimum becomes more broad, compared to the 
plots in figure 13). The pixel ages are consistent, however, with a picture where jet-induced 
star formation propagates outward from a region near the central "bar" in MO. 

We also extracted SEDs from the registered, convolved, binned multi-wavelength images 
in three boxes at the locations shown in Fig. 7, corresponding to MO-South, MO-North, and 
the filamentary structure to the north of MO. Following the same procedure as for the pixel 
age map, but using these larger boxes to improve signal-to-noise, we obtain ages (with 68 % 
confidence limits) of G.GG^o'g^ Myr, 7.86]'^q'2i Myr, and 7.00^033 Myr respectively for the three 
components. These more sensitive measurements support the conclusion from the pixel age 
map that MO-North was the first region to form stars, and that star formation spread over 
~ 1 Myr to the regions further out (MO South and the NE Filament have the same age 
to within these quoted errors). Presumably the star formation will continue to advance 
"downstream" in future, as the gas seen in the H I map begins to form stars. 



4. Discussion 

Our new observations of MO reinforce the conclusion by van Breugel et al. (1985) that 
MO is a young starburst triggered by the radio jet from NGC 541. However, some interesting 
questions remain - firstly, what existed at that location before the jet entered? Possibilities 
range from a pre-existing gas-rich galaxy, a cold gas cloud, or else a warm or even hot diffuse 
cloud of gas. Secondly, what are the physical conditions in the jet (speed, Mach number, 
mass loading etc.) such that it is able to trigger this starburst? Our observations add some 
further insight into the first question. To learn more about the second will probably first 
require a better understanding of FR-I jets in general (Young et al. 2005; Eilek et al. 2002; 
Kataoka et al. 2006) and, ultimately, numerical simulations of these types of jets. 
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4.1. Kinematics and morphology 

The neutral hydrogen emission is intimately associated with the radio jet: the H I 
appears to be "wrapped around" the axis of the jet, and there is no H I outside the jet! 
It will be important to search for possible very low surface brightness H I that might exist 
near MO, in the bridge or elsewhere, since our observations were taken in C-array and very 
large scale emission might have been missed. But given the current observational evidence, 
there is good reason to believe that there was no extended H I prior to the jet interaction, 
and therefore that the H I is entirely the result of radiative cooling of warmer gas in the 
inter-galactic medium triggered by the shocks induced in it by the passage of the radio jet. 

Numerical simulations of jet / cloud collisions, subject to assumptions about the radio 
jet and the multi-phase medium with which it interacts, show that radiative cooling may 
indeed lead to star formation (Fragile et al. 2004). The mass inferred for the H I cloud is 
4.9 X 10® Mq, and the projected extent of the H I emission is 275 kpc^ (section 2.7). The total 
amount of (hot) gas in this area according to the observations of Nikogossyan et al. (1999) 
is expected to be ~ 1.4 x 10^ Mq, so the available reservoir is certainly large enough. The 
question is: how did it cool to form the H I? Until we have better X-ray data we will have 
to assume that the IGM in the Abell 194 cluster resembles to some extent that of so-called 
cooling flow clusters, which contain both hot (10^ — 10^ K) and cool / warm (10'^ K) gas. We 
note Abell 194 is reportedly a "poor and cold" cluster (Nikogossyan et al. 1999) and as such 
a population of relatively cool, dense clouds in its IGM may well exist. 

The jet-H I relationship in MO is different from that in Centaurus A, where the jet is 
running into a pre-existing, very extended H I structure found in rotation around NGC5128 
(Schiminovich et al. 1994). Another difference is that in MO the H I is downstream and the 
star formation is closer to the AGN, while in Gen A it is the other way around: the H I 
is closest to the AGN and the star formation is downstream. Perhaps this is the difference 
between a jet hitting an existing H I cloud (Gen A) and the jet causing the warm gas to cool 
to form HI (MO). 

The H I shows a significant velocity gradient (up to 100 kms~^) at the jet / cloud 
interaction site (Oosterloo & Morganti 2005). Other examples, of even larger H I outflows 
(1000 kms~^) have recently been discovered in a number of radio galaxies, to the extent 
that these may remove significant amounts of gas from their parent galaxies (Morganti et 
al. 2005). However, these are all very luminous far-infrared galaxies with lots of ISM cold 
gas and dust through which the jets are propagating. In these cases it is reasonable to 
assume that the H I gas may have been entrained in some way, as it may also have in Gen A 
(Oosterloo & Morganti 2005). 
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In MO the velocity gradient is much smaller, ~ 40 kms~^, but the scenario may be 
qualitatively the same although in MO the entrained gas clouds were probably not cold but 
warm. However, also in this case the clouds will be accelerated, while they cool radiatively. 
Following Fragile et al. (2004) the velocity v^ of the accelerated clouds is 



kms^V ~ V103/ VlOOpc/ VIO^ kms-V ^Myr 



where x = f^ci/nb^i is the density contrast between the cloud and background intercloud 
medium; R^i is the cloud radius; Vsh,b is the shock velocity in the background intercloud 
medium, and t is the time during which the cloud is being accelerated. This time will 
be at most tacc ~ \fx^cc yr, where tec is the cloud crossing time for the radiative shocks. 
For the parameters used by Fragile et al. one finds tec ~ 1 Myr, tacc < 30 Myr and Vd ~ 
10 — 100 kms~^. Thus within the assumed framework of a jet interacting with a clumpy 
warm medium it is possible to explain the observed H I velocity shear of ~ 40 kms~^. 

The VLA H I, GALEX UV, and ATT Ha images all reveal that MO has a double 
structure which straddles the radio jet. In fact, the overall UV appearance of MO now 
looks very much like the radio-aligned rest-frame UV structure seen in the z = 3.8 radio 
galaxy 4C 41.17 (van Breugel et al. 1999), which is also thought to be due to jet-induced 
star formation (Dey et al. 1997; Bicknell et al. 2000). This an important qualitative result 
because it seems now more appropriate to assume that the star formation in MO occurs in 
the shear layers at the wall of the jet, i. e. not in a head-on jet-cloud collision. This means 
that the jet / cloud collision in MO may be a more complex astrophysical problem than 
initially thought. 

The jet speed in such a shear layer will be lower than at the center, due to jet expansion 
and entrainment (Laing & Bridle 2002). The pressure in the wall shock is also lower than 
in the bow-shock of the jet. This translates to a lower speed for the radiative shocks in the 
entrained clouds, and the parameter space studied by Fragile et al. (2004) then indicates 
that star formation may be easier to accomplish, at least with a certain range of pre-shock 
cloud densities. But perhaps more important is that shear layers of FR-I type jets are active 
sites of relativistic particle acceleration as evinced by radio spectral tomographic imaging 
which shows a remarkably fiat radio spectral energy distribution across the jets (Young et 
al. 2005). This is further supported by recent X-ray observations of the nearest FR-I radio 
jet, in Centaurus A (Kataoka et al. 2006), which show that the jet is edge-brightened with 
a possible hardening of the X-ray spectrum in the outer regions. As a result the pressure of 
the relativistic fiuid in these regions will be enhanced, which may affect the interaction with 
the clouds and the star formation process. 
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The peak H I column densities in the H I blobs are about Nhi = 4.5 x 10^*^ cmT^ . This 
is just at the bottom of the range where, under galactic conditions, star formation would 
occur (Schaye 2004), but considering the relatively low spatial resolution of the H I data, it 
could match more or less. However, the observed star formation sites are several kpc away 
from where the peaks in the H I data are, i. e. , the star formation is not happening where 
the H I says it should happen. So it may not be obvious to make the usual link between star 
formation and H I column density. The same is true in the H I filament near Centaurus A 
where the H I column densities are well above 10^^ cm~^ but only at those locations where 
there is no star formation. 

One can also compare the global star formation efficiency Mstars/Mui in MO to that 
of other systems. For MO we find ~ 4%, which is similar to that in Centaurus A (a few 
percent; Oosterloo & Morganti 2005), and indeed typical for starbursts in isolated dwarf 
galaxies (Fritze-v. Alvensleben 2005, and references therein). 

It is puzzling that, while the H I emission downstream from MO-South is comparable to 
that of MO-North, the Ha and UV emission and hence the current star formation activity 
is much less there. It should be noted in this regard that the jet is disrupted in the region 
of MO, so it may well be that the jet-cloud interaction is simply stronger in MO-North than 
in MO-South. This would be the case if, for example, the sense of the jet precession was 
towards the north. 

The HST images for the first time reveal the presence of numerous "blobs" which appear 
to be young clusters associated with H II regions. We extracted photometry in 0.5"-radius 
apertures for 20 of these blobs (Fig. 15), converting from the HST filters to Johnson V and 
/, and transforming the I^-band measurements to absolute magnitudes. Comparison with 
van den Bergh (1995) shows that the absolute magnitudes are what would be expected of 
young star clusters, and the colors of these blobs, when compared to Fig. 11 of Whitmore 
& Schweizer (1995), are just what would be expected for ages of a few Myr (consistent with 
all our other data). f/-band data would be necessary in order to pin down these ages more 
accurately, and we would need to reach fainter magnitudes in order to compare with the 
color-magnitude diagrams for individual stars in Centaurus A by Mould et al. (2000). 
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Fig. 15. — Left: F555W image of MO, showing the apertures used for F555W/F814W 
photometry of "blobs", numbered sequentially from south to north. Right: Photometry of 
the "blobs" in MO, converted from F555W/F814W to Johnson V and /, and to absolute 
magnitudes in V. Colors and magnitudes are typical for young star clusters. 
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4.2. Timing 

Due to the tendency of a young population to dominate tlie integrated starliglit in a 
galaxy, it is difficult to completely rule out the presence of an older pre-existing population. 
The situation in MO is even more complex due to its location in the outskirts of the elliptical 
galaxy NGC541, and in the bridge of stars connecting NGC541 with two of the other 
ellipticals in the cluster. Certainly, the spectroscopic and photometric data presented here 
agree that a population of order 7.5 Myr old is present in MO, and dominates the light output. 
The consistency of the pixel ages, the low internal reddening, and the lack of absorption lines 
expected from an old population, also support the conjecture that MO consists of a single-age 
population. Additionally, the similar morphologies at F555W and K' are not what would be 
expected if an old, gas-rich dwarf galaxy wandered into the jet and had a new burst of star 
formation "switched on" by the jet-cloud interaction. The star formation also appears at 
the leading edge of the H I cloud, as would be expected if it was the jet that caused the H I 
and subsequently the stars to form (since it takes time for the stars to form after the passage 
of the jet). As discussed above, however, it is not possible using these data to completely 
rule out the presence of an older underlying population in MO. It does seem clear, however, 
that the picture of a jet-cloud interaction, sometime a little more than 7.5 Myr ago, would 
be capable of causing the observed phenomena. 

Examination of the radio maps shows that a jet plume (van Breugel et al. 1985) extends 
up to 37kpc "downstream" of MO (or somewhat further if the jet is not in the plane of the 
sky). It seems unlikely that this is a buoyant bubble (Briiggen & Kaiser 2001) as the jet 
is still "switched on" and the plume does not display the characteristic "mushroom cloud" 
shape. The fact that the jet is well-collimated until its interaction with MO, and disrupted 
downstream, as well as the similar size of the H I cloud to the width of the jet, suggest 
that as the jet expanded, it ran into the gas cloud that ultimately formed MO. If MO was 
a pre-existing galaxy that wandered into the jet, it is unlikely that it would have done so at 
the precise time to cause the decollimation seen. 

The age of the downstream cocoon can therefore be estimated as the age of MO, plus 
the time taken for stars to begin forming after the jet-cloud interaction. Taking a typical 
star formation timescale < 1 Myr (Fragile et al. 2004), we assume that the cocoon expanded 
by at least 37kpc in 8.5 Myr, i. e. , at > 4 x 10^ kms"""^. 

If we assume that the SFR has been roughly constant during the lifetime of MO, this 
yields a timescale for the formation of the 1.9 x 10^ Mq of stars of 36.5 Myr. This is only a 
few times the 7.5 Myr inferred from the SED age dating. Since it is plausible that the SFR 
was somewhat higher in the past, it is reasonable to assume that all the stars may have been 
formed since the jet-cloud interaction. The timescale for depletion of the observed H I is 
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much longer (~ 1 Gyr). 



5. Conclusions 

We have presented new evidence that the star formation in MO was induced by a radio 
jet, primarily in the outer shear layers. This is similar to the jet-induced star formation 
associated with the Centaurus A jet, and the radio-aligned star forming regions in powerful 
radio galaxies at high redshift. Absorption and emission line measurements, and broad- 
band SED fitting, give a consistent age of around 7.5 Myr for MO. While it is not possible 
to completely rule out the presence of an old population in MO, the data are consistent 
with MO having formed de novo when the jet interacted with the ambient ISM / IGM. 
We discovered a double neutral hydrogen cloud which straddles the radio jet downstream 
of MO, at a location where the jet changes direction and de-collimates. Unlike Centaurus 
A, we propose that the H I associated with MO formed in situ through cooling of clumpy 
warm gas, in the stellar bridge or cluster IGM, as it was compressed by radiative shocks at 
the jet collision site, in agreement with earlier numerical simulations (Fragile et al. 2004). 
The star formation in MO then followed from the cooling and collapse of such H I clouds, 
and the H I kinematics, which show 40 kms~^ shear velocites, are also consistent with such 
models. The FR-I type radio source associated with MO is orders of magnitude lower in 
radio luminosity, and correspondingly more common, than the very luminous FR-II type 
radio galaxies and their aligned star forming regions at high redshift. Nevertheless, as MO 
shows, also these types of radio sources can trigger star formation, even if the ambient 
medium is only moderately dense. This suggests that jet-induced star formation may be 
quite common in the early universe, where gas densities in their parent, forming galaxies 
and surrounding proto-clusters were much higher than today, and when AGN activity was 
much more prevalent. 
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